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Abstract 

In this paper, we study the resource allocation algorithm design for distributed antenna multiuser 
networks with full-duplex (FD) radio base stations (BSs) which enable simultaneous uplink and downlink 
communications. The considered resource allocation algorithm design is formulated as an optimiza¬ 
tion problem taking into account the antenna circuit power consumption of the BSs and the quality 
of service (QoS) requirements of both uplink and downlink users. We minimize the total network 
power consumption by jointly optimizing the downlink beamformer, the uplink transmit power, and 
the antenna selection. To overcome the intractability of the resulting problem, we reformulate it as 
an optimization problem with decoupled binary selection variables and non-convex constraints. The 
reformulated problem facilitates the design of an iterative resource allocation algorithm which obtains 
an optimal solution based on the generalized Bender’s decomposition (GBD) and serves as a benchmark 
scheme. Furthermore, to strike a balance between computational complexity and system performance, a 
suboptimal algorithm with polynomial time complexity is proposed. Simulation results illustrate that the 
proposed GBD based iterative algorithm converges to the global optimal solution and the suboptimal 
algorithm achieves a close-to-optimal performance. Our results also demonstrate the trade-off between 
power efficiency and the number of active transmit antennas when the circuit power consumption is 
taken into account. In particular, activating an exceedingly large number of antennas may not be a power 
efficient solution for reducing the total system power consumption. In addition, our results reveal that 
FD systems facilitate significant power savings compared to traditional half-duplex systems, despite the 
non-negligible self-interference. 
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I. Introduction 

The next generation wireless eommunication systems are required to support ubiquitous and 
high data rate communieation applications with guaranteed quality of service (QoS). These re¬ 
quirements translate into a tremendous demand for bandwidth and energy consumption. Multiple- 
input multiple-output (MIMO) is a viable solution for addressing these issues as it provides extra 
degrees of freedom in the spatial domain which facilitates a trade-off between multiplexing gain 
and diversity gain. Hence, a large amount of work has been devoted to MIMO communication 
over the past decades (H, [|2l. However, the modest computational capabilities of mobile devices 
limit the MIMO gains that can be achieved in practice. An attractive alternative for realizing 
the performance gains offered by multiple antennas is multiuser MIMO, where a multiple- 
antenna transmitter serves multiple single-antenna receivers simultaneously Q, (H. In fact, the 
combination of multiuser MIMO and distributed antennas is widely recognized as a promising 
technology for mitigating interference and extending service coverage [15]|-[|7]|. Specifically, dis¬ 
tributed antennas introduce additional capabilities for combating both path loss and shadowing by 
shortening the distances between the transmitters and the receivers. Nevertheless, if the number 
of antennas is very large, the circuit power consumption of distributed antenna networks becomes 
non-negligible compared to the power consumed for transmission. However, this problem has 
not been considered in most of the existing literature [15]|-[|71 on power efficient communication 
network design. Furthermore, even with these powerful MIMO techniques, spectrum scarcity is 
still a major obstacle in providing high speed uplink and downlink communications. 

Traditional communication systems are designed for half-duplex (HD) transmission since 
this mode of operation facilitates low-complexity transceiver design. In particular, uplink and 
downlink communication are statically separated in either time or frequency, e.g. via time division 
duplex or frequency division duplex, which leads to a loss in spectral efficiency. Even though 
different approaches have been proposed for improving the spectral efficiency of HD systems, e.g. 
dynamic uplink-dowlink scheduling/allocation in time division duplex communication systems 
O, flU, the fundamental spectral efficiency loss induced by the HD constraint remains unsolved. 
On the contrary, full duplex (FD) transmission allows downlink and uplink transmission to occur 
simultaneously at the same frequency. In fact, FD radio has the potential to double the spectral 
efficiency of conventional HD communication systems. However, in practice, the downlink 
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transmission in FD systems creates self-interference to the uplink receive antennas which can be 
exceedingly large compared to the received power of the useful information signals. In fact, the 
huge difference in the power levels of the two signals saturates the dynamic range of the analog- 
to-digital converter (ADC) essentially preventing FD communication. Fortunately, several recent 
breakthroughs in hardware (/signal processing algorithm) design for suppressing self-interference 
have been reported and FD radio prototypes have been successfully presented [[T0]| - lfT4]l . As a 
result, FD radio has regained the attention of both industry ifTSll - lfT^ and academia [[T^ - [|2^ . 
In lfT9ll . the authors studied techniques for self-interference suppression and cancellation for FD 
multiple-antenna relays. In ll20l . the outage probability of MIMO FD single-user relaying systems 
was investigated. In [l2T]l . a resource allocation algorithm was proposed for maximization of the 
achievable end-to-end system data rate of multicarrier multiuser MIMO FD relaying systems. 
In [l22l . a suboptimal beamformer design was considered to improve the spectral efficiency of a 
FD radio base station enabling simultaneous uplink and downlink communication. In lf2^ . the 
concept of FD communication was extended to the case of massive MEMO where a FD radio 
relay is equipped with a large number of antennas for suppressing the self-interference and for 
enhancing the system throughput. However, the benefits of multiple-antenna FD radio do not 
come for free. The rapidly escalating cost caused by the power consumption of the circuitries of 
large antenna systems has lead to significant financial implications for service providers, which is 
often overlooked in the literature [fTOl - ll^ . In fact, the systems in [[T0ll - lf23]l are designed to serve 
peak service demands by activating all available antennas of the system, without considering the 
power consumption in the off-peak periods. However, the service loads vary across a wireless 
network in practice, depending on the geographic location of the receivers and the time of day. 
Thus, we expect that extra power savings can be achieved by dynamically switching off some of 
the antennas. Nevertheless, the optimal number of active antennas has not been investigated from 
a system power efficiency point of view for FD radio communication, yet. In addition, there 
may be fewer degrees of freedom for self-interference suppression at each FD radio base station 
in distributed antenna systems if the total number of antennas in the network is fixed. Thus, 
it is unclear whether the distributed antenna architecture leads to power savings for FD radio 
communication. Furthermore, unlike for the orthogonal transmission adopted in HD systems, 
the uplink and downlink transmit powers are coupled in FD systems which make the design of 
efficient resource allocation algorithms particularly challenging. 


4 


In this paper, we address the above issues and study the resouree alloeation algorithm design for 
multiuser distributed antenna eommunieation networks. We minimize the total network power 
eonsumption while taking into aeeount the eireuit power eonsumption of the distributed BS 
antennas and ensuring the QoS of both uplink and downlink users. In partieular, we propose an 
optimal iterative resouree alloeation algorithm based on the generalized Bender’s deeomposition 
m-m- Furthermore, we propose a suboptimal resouree alloeation seheme with polynomial 
time eomputational eomplexity based on the differenee of eonvex funetions (d.e.) programming 
[l27]l whieh finds a loeal optimal solution for the eonsidered optimization problem. 

II. System Model 

A. Notation 

Matriees and veetors are represented by boldfaee eapital and lower ease letters, respeetively. 
A^, A^, Tr(A), and Rank(A) represent the Hermitian transpose, the transpose, the traee, and 
the rank of matrix A, respeetively; A >- 0 and A ^ 0 indieate that A is a positive definite and a 
positive semidefinite matrix, respeetively; Ijv is the NxN identity matrix; and denote 

the sets of all A x M matriees and NxN Hermitian matriees with eomplex entries, respeetively; 
diag(a;i, • • • , xk) denotes a diagonal matrix with the diagonal elements given by {xi, • • • , xk}', 
I'l denotes the absolute value of a eomplex sealar; the eireularly symmetrie eomplex Gaussian 
distribution is denoted by CAA(^, C) with mean veetor and eo-varianee matrix C; ~ stands for 
“distributed as”; £{■} denotes statistieal expeetation; and Vx/ denotes the gradient of a funetion 
/ with respeet to veetor x. 


B. System Model 

We eonsider a distributed antenna multiuser eommunieation network. The system eonsists of 
a eentral proeessor (CP), L FD radio base stations (BSs), and K mobile users, ef. Figure [IJ Eaeh 
FD radio BS is equipped with At > 1 antennas for downlink transmission and uplink reeeptiorH. 
The K users employ single-antenna HD mobile eommunieation deviees to ensure low hardware 
eomplexity. In partieular, and K-q users are seheduled for simultaneous uplink and downlink 
transmission, respeetively, sueh that Au + Ad = A. On the other hand, the CP is the eore unit of 
the network. In partieular, the FD radios are eonneeted to the CP via baekhaul links. In addition. 


'We assume that the antennas equipped at the FD BSs can transmit and receive simultaneously which has been successfully 
demonstrated in some FD radio prototypes fT2l . 





5 


FD radio BS 3 


Signal 


Co-channel Interference 


Central processor 
(CP) 


Backhaul connection 


Fig. 1. Multiuser downlink distributed antenna communication system model with L = 3 full duplex (FD) radio base stations 
(BSs), K\j = 1 uplink user, and Ku = 1 downlink user. For the depicted case, the antennas equipped at FD radio BS 2 are 
switched to idle mode for reducing the total power consumption in the network. 

the CP has the full channel state information of the entire network and the data of all downlink 
users for resource allocation. In this paper, we assume that the CP is a powerful computing unit, 
e.g. a series of baseband units as in cloud radio access networks (C-RAN), which computes the 
resource allocation policy and broadcasts it to all FD radio BSs. Each FD radio BS receives the 
control signals for resource allocation and the data of the downlink users from the CP via a 
backhaul link. Furthermore, the FD radio BSs transfer the received uplink signals via backhaul 
links to the CP, where the information is decoded. In this paper, we assume that the backhaul 
links are implemented with optical fiber and have sufficiently large capacity and low latency to 
support real time information exchange between the CP and the FD radio BSs. For studies on 
the impact of a limited backhaul capacity on the performance of wireless systems, please refer 

to m, iHHi. 

C. Channel Model 

A frequency flat fading channel is assumecj^ in this paper. The received signals at downlink 
user k G Kd} and the L FD radio BSs are given by 



and 


( 1 ) 


CO—channel interference 


y 
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self—interference 


+Z, 


( 2 ) 


The frequency flat fading channel can be interpreted as one subcarrier of an orthogonal frequency division multiplexing 


system. 
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respectively, where x G denotes the joint transmit signal vector of the L FD radio BSs to 

the Kj) downlink users. The downlink channel between the L FD radio BSs and user k is denoted 
by hoj. G ^ G C to represent the channel between uplink user j and 

downlink user k. and are the transmit data and transmit power sent from uplink user j to 
the L FD radio BSs, respectively, hu^. G uplink channel between uplink user j and 

the L FD radio BSs. Due to simultaneous uplink reception and downlink transmission at the FD 
radio BSs, self-interference from the downlink impairs the uplink signal reception. In practice, 
different interference mitigation techniques such as antenna cancellation, balun cancellation, and 
circulators ifT^ . [fTSll have been proposed to alleviate the impairment caused by self-interference. 
In order to isolate the resource allocation algorithm design from the specific implementation of 
self-interference mitigation, we model the residual self-interference after interference cancellation 
by matrix Hsi G Variables Iid^,, gj,k, Hsi, and hu^ capture the joint effect of path 

loss and multipath fading, z ~ CAA(0, ctIIattl) and rik ~ CAA(0, represent the additive white 
Gaussian noise (AWGN) at the L FD radio BSs and user k, respectively. 

In each scheduling time slot, Kb independent signal streams are transmitted simultaneously at 
the same frequency to the Kb downlink users. Specifically, a dedicated downlink beamforming 
weight, wl G C, is allocated to downlink user k at the l-th, I G {1, ..., antenna to 

facilitate downlink information transmission. For the sake of presentation, we define a super¬ 
vector Wfc G for downlink user k as 

= [wlwl ... . ( 3 ) 

Wfc represents the joint beamformer used by the N^L antennas shared by the FD radio BSs for 
serving downlink user k. Then, the information signal to downlink user k, x^, can be expressed 
as 

Xfc = (4) 

where G C is the data symbol for downlink user k. Without loss of generality, we assume 
that S{\dm = = 1, Vfc G {1,..., Kb},j e {1,..., Kb}. 


D. Network Power Consumption Model 

In our system model, we include the circuit power consumption of the system in the objec¬ 
tive function in order to design a resource allocation algorithm which facilities power-efficient 
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communication. Thus, we model the power dissipation in the system as the sum of one static 
term and four dynamie terms as follows lf2^ : 

N'j' L L 

ZYtp (w,, Si, = Po + ^ 

1=1 1=1 

" -V-' 

Antenna power eonsumption 

NtL Xd Ku 

+ ^^uCjP/) (5) 

1=1 k=l j=l 

" -V-' 

Amplifer power eonsumption 

where Pq is the aggregated static power eonsumption of the CP, all FD radio BSs, and all baekhaul 
links. Si E {0, 1} is a binary seleetion variable. In partieular, s; = 1 and s; = 0 indicate that 
the /-th antenna in the FD eommunieation system is in aetive mode and idle mode, respeetively. 
Si will be optimized to minimize the total network power eonsumption in the next seetion. 
pActive > 0 is the signal proeessing power that is eonsumed if an antenna is aetive. 
ineludes the power dissipations of the transmit filter, mixer, frequeney synthesizer, digital-to- 
analog converter, etc. In this paper, an FD radio antenna is eonsidered aetive if it serves at least 
one user in the system, pi'll® > 0 is the required power consumption of an antenna in idle mode, 
i.e., if it is not serving any user, and > pwie holds in general. Yl!k=i Yld=i\'^k\^ *^he 

total power radiated by the L FD radio BSs for downlink transmission. £d > 1 and Su > 1 
are eonstants whieh aeeount for the ineffieieney of the power amplifier adopted for downlink 
and uplink transmission, respeetively. In other words, Sd Yl!k=i Yl!i=i\'^k\^ Yl!j=i Pf 

the total power eonsumptions of the power amplifiers for downlink and uplink transmission, 
respeetively. > 0 and Q > 0 in the last two terms of dS]) are constant weights which can be 
ehosen by the system designer to prioritize the importanee of the total downlink transmit power 
and the transmit power of individual uplink users j E iCu}> respectively. 

III. Problem Formulation 

In this seetion, we first introduee the QoS metries for the considered FD radio eommuni- 
eation network. Then, we formulate the resouree allocation algorithm design as a non-eonvex 
optimization problem. 


^We assume Class A power amplifiers with linear characteristic are implemented in the transceivers. In practice, the maximum 
power efficiency of Class A amplifiers is 25%. 






A. Achievable Data Rate 


The achievable data rate (bit/s/Hz) between the L FD radio BSs and downlink user k G 
{1,... ,iTD} is given by 


Ck = log 2 (l + Tj^ ), where = 


Kn 


( 6 ) 




t^k 


2 

life 


is the receive signal-to-interference-plus-noise ratio (SINK) at downlink user k. 

On the other hand, we assume that the CP employs a linear receiver for decoding of the 
received uplink information. Therefore, the achievable data rate between the L FD radio BSs 
and uplink user j is given by 


CE = iog2(i + r 


UL 

3 


Kn 


rUL 

(7) 

Kv 


Hg)v,+y:/-"|vfhu.r. 

(8) 


k=l r^j 

where Vj G is the receive beamforming vector for decoding of the information for uplink 

user j. In this paper, maximum ratio combining (MRC) is adopted, i.e., the receive beamformer 
for uplink user j is chosen as = E^^i^ SiRihuj to maximize the signal strength of the received 
signal, where R; = diag ( 0, • • • , 0,1, 0, • • ■ , 0 ), V/ G {1,..., LNt}, is a diagonal matrix. It is 


(l — l) LN'y — I 

known that MRC achieves a good system performance, especially if a large number of antennas 
is employed, and has been widely adopted in the literature ll^. ||4l. 

Remark 1: We note that zero-forcing beamforming (ZFBF) or minimum mean square error 
beamforming (MMSE-BF) are not considered for uplink signal detection since they do not 
facilitate an efficient resource allocation algorithm design for the considered network. 

Using MRC, the uplink SINR of user j is given by 


pUL ^ 
^ 3 




Pf Tr (^h„,hg. Y.7.\ SSf R: 


H 


where 


Ky) N^L 

Tr ( ^ HsiWfcwf ^ ^ s^s„R„hu^R(f.R: 

k=l m=l n=l 

N’y L A^x ^ 

+ Pr Tr (hu,h{f^ SmSnRmhu.h^^R 

r^j m=l n=l 


H 


(9) 

( 10 ) 

( 11 ) 





B. Optimization Problem Formulation 


The system objeetive is to minimize the total network power consumption while providing 
QoS for reliable communication to both uplink and downlink users simultaneously. We obtain 
the optimal resource allocation algorithm policy by solving the following optimization problem: 


s.t. Cl: 


minimize Wtp w^, 

V ^ 


Ku 


>r? 4 , vfcG{i,...,PD}, 


E Pf\93,k\^ + 

t^k 

C2: ryL>rUL Vje iCu}, 


2 

life 


J 

Kt> 


C3: ^ Oi-Pmk.'" e {1. ■ ■ ■. IVxi}, C4: 0 < < Fl,., Vj e {1,..., A'„}, 

( 12 ) 


k=l 


C5: 5ze{0,l},V/e ATtL}. 


and ryy in constraints Cl and C2 denote the minimum receive SINK required by downlink 

^AC 

user k and uplink user j for successful information decoding, respectively. In C3, we constrain 
the maximum radiated power of the /-th antenna in the system to satisfy the maximum 

power spectral mask limit. C4 limits the maximum transmit power and ensures the non-negativity 
of the transmit power of uplink user j. C5 constrains the optimization variables which control 
the active and idle states of the antennas in the system to be binary. 

Remark 2: In this paper, energy/power saving is achieved by optimizing not only the uplink 
and downlink transmit powers, but also by optimizing the states of the antennas in the network. 
Thereby, it is expected that switching the antennas on and off adaptively according to the channel 
conditions is an effective strategy for reducing the network power consumption when the QoS 
requirements are not stringent or the number of users is low. 


IV. Resource Allocation Algorithm Design 

The optimization problem in (O is a mixed non-convex and combinatorial optimization 
problem. The combinatorial nature is due to the binary selection variables in C5. Also, variable 
Si is coupled with both downlink beamforming vector and uplink power allocation variable 
in constraint C2. Furthermore, constraint Cl is non-convex with respect to w^. In the 
following, we first transform the optimization problem into an equivalent form and obtain the 
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global optimal solution by using the generalized Bender’s decomposition. Then, we propose a 
suboptimal polynomial time algorithm which is inspired by the difference of convex functions 
program. 


A. Problem Reformulation 

In this section, we reformulate the considered optimization problem in (fT^ using the defini¬ 
tions Wfc = Wfcwf, and Hu^. = hu^.h^,. This leads to 


Kv 


N^L N'-^L K-q 

minimize Pq + ^ + ^(1 - + r/ Tr(Wfc) + eu Ci-P, 


V 


i=i 


Cl: 3hh£^>^Tr(Ho,Wj) + |^p«|g,.tP + <,Vi, 

t^k j = l 

p\] A^t-^ N'yL 

C2: Tr THu^. ^ ^ gm,nR-m,Huj.R^ 

^ reqj 

N^L A^x-^ A^x^ 

> Tr ( ^ gMHugRz) + Tr ( HsiW^Hg ^ ^ g^,„R„Hu,R^ 

1=1 k=l m=l n=l 

A^xA A^xA 

+ Tr (Hu, ^ ^ pUg,^,„R„Hu,R^), Vj G {1,..., K^}, 

r^j m=l n=l 

A'd 

C3: 5^Tr(W,R0<szP°L^,V/e {!,..., iVTL}, C4, C5, 

fc=i 

C6: Wfc P 0, Vfc, C7: Rank(Wfc) < 1, V/c, C8: 0 < qm,n < Sm, Vm, n G {1,..., A^tP}, 


Qm^n — ^ni VlTT, Tl^ CIO. qrn,n P '^n “f f: '^TTl^ Tl. 


(13) 


Constraints C6, C7, and G are imposed to guarantee that = Wkwjl holds 

after optimization, g^.n is an auxiliary continuous optimization variable which is introduced to 
handle the product of two binary variables SnSm in constraint C2, cf. (l^- (fTT]) . In particular, 
because of constraints C8 - CIO, qm,n will have a binary value if si is binary. 

We note that constraint C2 is still non-convex due to the product terms qm,nP^ and qm,ny^k 
which is an obstacle for the design of a computationally efficient resource allocation algorithm. 
In order to circumvent this difficulty, we adopt the big-M formulation ll29l . [[30ll to decompose 
the product terms. First, we introduce auxiliary variables = Pj^<lm,n and W™’” = Wkqm,n- 
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Then, we impose the following additional constraints: 


Cll: P. 


.u 

j,m,n 




C12: PL„ <pu (14a) 


C13- P’^ > P'^ 

j,m,n — ^ 3 


(1 - gm,n)Pmax,-, Vj, Trt, n, C14: > 0, 


3,m,n 

u 


(14b) 


C15: ^ WP^x^W, VA:,m,n, 


C16: W™’” P Wfc, VA:, m, n, (14c) 


C17: W™’" P Wfc - (1 -g^,„)lAr^iP°L vfc,m,n,C18: W™’" P 0,Vfc,m,n. 


(14d) 


In particular, constraints C11-C18 involve only continuous optimization variables, i.e., Pj, Pj,m,n, 
qm,n, and W^, which facilitates the design of an efficient resource allocation algorithm. Subse¬ 
quently, we substitute P^^^.n = Pfqm,n and W™’” = Wkqm,n into the coupled variables in C2 
which yields 


C2: 


> 


rUL 

req^. 


A^X-^ A^X-^ 


Tr (Hu, Y. E ^E,..H.»Hu,R 


m=l n=l 

N'yL A^d A^xA N'yL 

uj 'll' ( E «mH„,R,) + tv (E E E HsiW”’"H«R„H„,R: 

/=1 k=l m=l n=l 


H 


N'j^L A^xA 


(15) 


+ E Ti' (Hu, E E P»..»R».Hu^R"j, Vj e {1,..., Au}. 

r^j m.=l n=l 

The big-M formulation linearizes the terms qm,nP^ and qm,ny^k such that constraint C2 is an 
affine function with respect to the new optimization variables P^rn,n and W^’”. We note that 
constraints C2 and C2 are equivalent when constraints C5 and C11-C18 are satisfied. 

As a result, the considered optimization problem (fTSl) can be transformed into the following 
equivalent problem: 

A^xA N^L Ky) 

mininiize Po + V + V(1 - s^)P^^^" + V £d Tr(W;t 


pU pU 

j ’ j,m,n 


1 = 1 


1 = 1 


k=l 


Kv 

E'^oCi-p,' 

i=i 


s.t. 


C1,C2,C3, C4, C6,C8 - C18, 

C5: Si e {0,1}, V/, C7: Rank(Wfc) < 1, VA:, 


(16) 


and we can focus on the design of an algorithm for solving the optimization problem in (fT^ . 
Now, the remaining non-convexity of optimization problem (fT^ is due to constraints C5 and 
Cl. 

Remark 3: We note that the uplink-downlink duality approach in llSTI . Il3^ cannot be applied 
to our problem for the following two reasons. First, the uplink and downlink transmit power 
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variables are coupled in constraints Cl and C2. Second, the uplink and downlink transmit powers 
of each transceiver are constrained. 

B. Optimal Iterative Resource Allocation Algorithm 

Now, we adopt the generalized Bender’s decomposition (GBD) to handle the constraints 
involving binary optimization variables It2^ - lt26l . i.e., C3, C8, C9, and CIO. In particular, we de¬ 
compose the problem in (fT^ into two sub-problems: (a) a primal problem which is a non-convex 
optimization problem involving continuous optimization variables {W^, W^’", Pf, Ppm,n^ 

(b) a master problem which is a mixed integer linear program (MILP). Specifically, the primal 
problem is solved for given si which yields an upper bound for the optimal value of (fT^ . In 
contrast, the solution of the master problem provides a lower bound for the optimal value of (fT^ . 
Subsequently, we solve the primal and master problems iteratively until the solutions converge. 
In the following, we first propose algorithms for solving the primal and master problems in the 
i-th iteration, respectively. Then, we describe the iterative procedure between the master problem 
and the primal problem. 

1) Solution of the primal problem in the i-th iteration: For given and fixed input parameters 
Si = si{i) obtained from the master problem in the z-th iteration, we minimize the objective 
function with respect to variables {W^, W™’”, in the primal problem: 

NtL NtL Ko Kv 

minimize Pq + Y" + V(1 - + r? V £d Tr(W^) + £u V CjP^ 

Wj.gH t ,W^ 

p\J pU 
j ’ j,m,n 

s.t. Cl,C2,C3, C4, C16 - C18. (17) 

We note that constraint C5 in (fT^ will be handled by the master problem since it involves only 
the binary optimization variable si. Now, the only obstacle in solving (fTTI) is the combinatorial 
rank constraint in C7 and we adopt the SDP relaxation approach to handle this non-convexity. In 
particular, we relax constraint C7: Rank(Wfc) < 1 by removing it from the problem formulation, 
such that the considered problem in (fTTI) becomes a convex SDP and can be solved efficiently 
by numerical methods designed for convex programming such as interior point methods [|3^ . 
If the solution of the relaxed version of (fTTI) is a rank-one matrix for all downlink users, 
then the problem in (fTTI) and its relaxed version share the same optimal solution and the same 
optimal objective value. 
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Now, we study the tightness of the adopted SDP relaxation. The SDP relaxed version of 
(fTTl) is jointly eonvex with respeet to the optimization variables and satisfies Slater’s eonstraint 
qualifieation. Thus, strong duality holds and solving the dual problem is equivalent to solving 
(fTTl) . To obtain the dual problem, we define the Lagrangian of the relaxed version of (fTTl) as 


£(©,$) =WTp(Wfc,sz,pU)+/i(0,$) + / 2 (©,$), where 


3 

NtL 


( 18 ) 


N^L -^u 

ZYtp(w,, s;, Pf) =Po + J2 + Tr(Wfe) +£u (19) 

1=1 1=1 k=l j=l 

Ku Ku NtLNt^L iTu 

/.(©,«) = - Ti'(z,w^) - ^ Xi(pf - Pi. 

j=l m=l n=l j=l 

Ku 

k=l 


?U 


k=l 

T1(Hd.W, 


pDL 

reqfc 


JVd Ku Ku 

■- £iv(HDWi) + E P - E xjp 

t^k j=l j=l 


Kr 


+ ‘^-Tr y]g,.,H„,R, 


rUL 

req,. 


i=i 

L A^x ^ 

( E E E HsiW™'”HgR„H„,R") + E(Hu. E E P»..»R™Hu,R 


1=1 

A^x^ A^x^ 


H 


k=l m=l n=l 
Kjj A^x ^ -^T ^ 


m=l n=l 
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D 


CiSi. 


A'd A^tA A^xA A^xA 

M @,«) = E E ft (ti(WiR,) - + E E 

k=l 1 = 1 772=1 72=1 

A^xA N'-^L A^xA A^xA 

“1“ ^ ^ ^ ^ V^772,72 (Q'772,72 '^Tl) H“ ^ ^ ^ ^772,72 ('^72 ^771 1 0^772,72)* 


72 ’^mj 


( 21 ) 


771=1 72=1 772=1 72=1 

Here, © -{AV k-} SI ■} Pj 5 j Pj,m,n'> Qm,n} 3.11(1 $ 5 pl 5 7 Xj } 7 {^ 772,727 ^ 772 , 72}7 Pm,n’) 

/^J,772,72"^7,772,725 772,725 (^j.,m,n'} Dci5 ; 19ci 6, ; i l^Cis } are the eolleetions of 

’ ' ! ,z> ? j) 1 ji 1 ■^'^k,m,n ^^k,m,71 ^‘k,m,n ■^°k,m,n^ 






14 


primal and dual variables, respectively; > 0,p/ > 0, {Aj, Xj} > 0, ^ 0,<^m,n > 

0,Dci6 ^ 0,Dci 7 ^ 0, and Dcis — 0, are the scalar/matrix dual variables for 

°fc,m,n k,m,n ^°k,7n,n 

constraints Cl - C4, C8 - C18, respectively. Function si, in (fT9l) is the objective 

function of the SDP relaxed version of problem (fTTI) : /i(©, $) in (l20l) is a function involving 
the constraints that do not depend on the binary optimization variables; / 2 (©,$) in (l2TI) is 
a function involving the constraints including These functions are introduced here for 
notational simplicity and will be exploited for facilitating the presentation of the solutions for 
both the primal problem and the master problem. 

For a given si, the dual problem of the SDP relaxed optimization problem in (fTTI) is given by 


maximize minimize £(©,$). (22) 

We define ©*(*) = {W*,sz,pu*, and $(*) = {$*} as the optimal primal 

solution and the optimal dual solution of the SDP relaxed problem in (fTTI) in the i-th iteration. 

Now, we introduce the following theorem regarding the tightness of the adopted SDP relax¬ 
ation. 

Theorem 1: Assuming the channel vectors of the downlink users, hn^., k E {1,..., Pd}, can 
be modeled as statistically independent random variables, then the solution of the SDP relaxed 
version of (fTTf) is rank-one, i.e., Rank(Wfc) = l,V/c, with probability one. Thus, the optimal 
downlink beamformer for user k, i.e., w^, is the principal eigenvector of W^. 

Proof: Please refer to Appendix A for a proof of Theorem 1. 

On the other hand, we formulate an /i-minimization problem for the case when (fTTf) is 
infeasible for given binary variables si{i). The /i-minimization problem is given as: 


NtL 




minimize 




m,n ' m,n ' m,n 


1=1 


m=l n=l 


S.t. 


C1,C2,C4, C6, Cll - C18, 


(23) 


C3: Tr(W,R0 < G {1,..., N^L}, 

C8: 0 < qm,n < Sm{i) + Vm, n e {1,... , At Tv}, C9: qm,n < Sn{i) + Vm, n, 
CIO: + qm,n > Sn{i) + Sm{i) - 1, Vm, n, C19: z/p, z/^^, z/^^, z/^^ > 0, V/, m, n, k. 
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Equation (|2^ is an SDP problem and can be solved by interior point methods with polynomial 
time computational complexity. We note that the objective function in (|2^ is the sum of the 
constraint violations with respect to the problem in (fT71) . Besides, the corresponding dual variables 


and the optimal primal variables will be used as the input to the master problem for the next iter¬ 
ation lf25]l . We adopt a similar notation as in (fT3l) to denote the primal and dual variables in (|2^ . 


In particular, the primal and dual solutions for the /i-minimization problem in (12^ are denoted as 

® ky Si, P j y Pj^rn,ni Snd $ 'Ip j, Pi, {\j, Xj} y '^ky ^m,ny f^m,ny'^rn,ny ^rn,ny 

Hm,ny , Dc,,, ^J, respectively. The primal 

and dual variables will be exploited as inputs for the constraints of the master problem. 


2) Solution of the master problem in the i-th iteration: For notational simplicity, we define P 


and X as the sets of all iteration indices at which the primal problem is feasible and infeasible, 
respectively. Then, we formulate the master problem which utilizes the solutions of (fT3]) and 
(1231) . The master problem in the i-th iteration is given as follows: 


minimize g 


(24a) 




s.t. C5, 


(24b) 

Q > e {1,.. 

. ,i} nP, 

(24c) 

0 > e(^(t),sz),f e 

.,*}nx, 

(24d) 


where si and q are optimization variables for the master problem and 

= minimize Wtp (w^, sz, $(t)) +/^(e, (25) 

pU pU 

3 ’ j,m,n 

^($(t),sz) = minimize /i(0, $(t))-f/ 2 ( 0 , $(f)). (26) 

pU pU 

j ' j,m,n 

Equations (1251) and (1261) are two different minimization problems defining the constraint set 
of the master problem in (12^ . In particular, g > ,^($(f), si),t e {1,... ,i} D P in (I24cl) and 
0 > ^($(f),sz),f G {l,...,i}nX in (124dl) denote the sets of hyperplanes spanned by the 
optimality cut and the feasibility cut from the first to the Pth iteration, respectively. The two 
different types of hyperplanes reduce the search region for the global optimal solution. Moreover, 
both ^(^>(t), s;) and ^($(t), s;) are also functions of s; which is the optimization variable of the 
outer minimization in (1241) . 






TABLE I 


Optimal Iterative Resource Allocation Algorithm based on GBD 


Algorithm Generalized Bender’s Deeomposition 

1: Initialize the maximum number of iterations Lmax, UB(0) = oo, LB(0) = —oo, and a small constant 0 
2: Set iteration index i = 1 and start with siii) = 1, Vfc, I 
3; repeat {Loop} 

4; Solve O by SDP relaxation for a given set of si{i) 

5: if i ll7b is feasible then 

6; Obtain an intermediate resource allocation policy &'{i) = , W}™’" , Pj^m,ny <lrn,n}^ th® corresponding 

Lagrange multiplier set ^(i), and an intermediate objective value /q 
7: The upper bound value is updated with UB(i) = min{UB(i — l),/o}- If UB(i) = /q, we set the current optimal 

policy ©current = 

8; else 

9; Solve the -minimization problem in j23b and obtain an intermediate resource allocation policy &{i) = 

Qrn.n} and the corresponding Lagrange multiplier set $(j) 

10: end if 

11: Solve the master problem in l l24b for si, save si{i -|- 1) = si, and obtain the i-th lower bound, i.e., LB(i) 

12: if |LB(i) - UB(i)| < 1 ? then 

13: Global optimal = true, return {Wl, sf , P ^*= {©current} 

14: else 

15: i = i -I- 1 

16: end if 

17. until i — Lmax 

Now, we introduce the following proposition for the solution of the two minimization problems 
in (1251) and (l26l) . 

Proposition 1: The solutions of (1251) and (|2^ for index t G {1,... ,i} are the solutions of 
(fTTI) and (|2^ in the f-th iteration, respectively. 

Proof: Please refer to Appendix B for a proof of Proposition [T1 
The master problem in (l24l) is transformed to an MILP by applying Proposition [T] to solve 
(1251) and (|2^ . Hence, the master problem can be solved by using standard numerical solvers for 
MILPs such as Mosek ll34ll and Gurobi lIMll . We note that an additional constraint is imposed 
to the master problem in each additional iteration, thus the objective value of (l24l) . i.e., the 
lower bound of (fT^ . is a monotonically non-decreasing function with respect to the number of 
iterations. 

3) Overall algorithm: The proposed iterative resource allocation algorithm is summarized 
in Table U and is implemented by a repeated loop. For the initiation, we first set the iteration 
index i to one and the binary variables sfi) to one, e.g. 5^(1) = 1,V/. In the i-th iteration, we 
solve the problem in (fTTI) via SDP relaxation. Two different types of Lagrange multipliers are 
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defined depending on the feasibility of the primal problem. If the problem is feasible for a given 
si{i) (lines 6, 7), then we obtain an intermediate resouree alloeation poliey ©(1), an intermediate 
objeetive value /q, and the eorresponding Lagrange multiplier set $(1). In partieular, $(*) is used 
to generate an optimality cut in the master problem. Also, the optimal resouree alloeation policy 
and the performance upper bound UB(7) are updated if the computed objective value is the lowest 
across all the iterations. On the contrary, if the primal problem is infeasible for a given si{i) 
(line 9), then we solve the /i-minimization problem in (|2^ and obtain an intermediate resource 
allocation policy ©(i) and the corresponding Lagrange multiplier set $(1). This information will 
be used to generate an infeasibility cut in the master problem. We note that the upper bound 
is obtained only from the feasible primal problem. Subsequently, we solve the master problem 
based on ©(f) and ©(*), t G {1,... ,1}, via a standard MILP numerical solver. Due to weak 
duality the optimal value of the original optimization problem in (fTTI) is bounded below 
by the objective value of the master problem in each iteration. The algorithm stops when the 
difference between the Lth lower bound and the 1-th upper bound is smaller than a predefined 
threshold > 0 (lines 12 - 14). We note that when the master and the primal problems can be 
solved in each iteration, the proposed algorithm is guaranteed to converge to the optimal solution 
[l25l Theorem 6.3.4]. 

C. Suboptimal Resource Allocation Algorithm Design 

The optimal iterative resource allocation algorithm proposed in the last section has a non¬ 
polynomial time computational complexity due to the MILP master problem 0. In this section, we 
propose a suboptimal resource allocation algorithm which has a polynomial time computational 
complexity. The starting point for the design of the proposed suboptimal resource allocation 
algorithm is the reformulated optimization problem in (fTBl) . 

1) Problem reformulation via difference of convex functions programming: The major obstacle 
in solving (flSI) are the binary constraints. Hence, we rewrite constraint C5 in its equivalent form: 

A^x L A^x L 

C5a: 0 < s; < 1, V( e {1..., L} and C5b: sf < 0. (27) 

i=i i=i 

Now, optimization variables si in C5a are continuous values between zero and one while 
constraint C5b is the difference of two convex functions. By using the SDP relaxation approach 

''The optimal algorithm serves mainly as a performance benchmark for the proposed suboptimal algorithm. 


TABLE II 


SuBOPTiMAL Iterative Resource Allocation Algorithm 
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Algorithm Successive Convex Approximation 

I: Initialize the maximum number of iterations Lmax, penalty factor (f) ^ 0, iteration index i = 0, and 
2; repeat {Loop} 

3; Solve d31b for a given and obtain the intermediate resource allocation policy {W}, 

4: Set = sj and i = i + 1 

5: until Convergence or i = Lmax 


as in the optimal resource alloeation algorithm, we ean reformulate the optimization problem as 

rninimize_ Utp fw^, si, pf) (28) 

s.t. & eV, C5b, 

where V denotes the eonvex feasible solution set spanned by eonstraints Cl, C2, C3, C4, C5a, C6, 
and C8 - Cl8. The only non-convexity in (|2^ is due to eonstraint C5b whieh is a reverse eonvex 
funetion ll27]l . Now, we introduee the following Theorem for handling the eonstraint. 

Theorem 2: For a large eonstant value (j) ^ 1, (l28l) is equivalenj^ to the following problem: 

minimize Utp (w^, sz, pf) + </> ( ^ 

S.t. & ev. 

In partieular, 0 aets as a large penalty faetor for penalizing the objeetive funetion for any si that 
is not equal to 0 or 1. 

Proof: Please refer to Appendix C for a proof of Theorem [2l 
The problem in (|2^ is in the canonieal form of differenee of eonvex (d.c.) funetions pro¬ 
gramming. Speeifieally, g[si) = sf is a eoneave funetion and we minimize d.e. funetions 

over a eonvex constraint set. As a result, we ean apply sueeessive eonvex approximation to 
obtain a local optimal solution of (l29l) . 

2} Suboptimal iterative algorithm: Since g{si) is a differentiable convex function, inequality 

9{si) > g{sP) + V,M4''){si-sf'>),VlE{f...,N^L}, (30) 

always holds for any feasible point sf\ where the right hand side of (l30l) is an affine funetion 
[l3^ and represents a global underestimator of g{si). 


5 


Here, equivalence means that both problems share the same optimal objective value and the same optimal resource allocation 


policy. 
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As a result, for any given value of we solve the following optimization problem, 


minimize 


s.t. 


N'j' L N'-£ L N'x L 

(0'\2 o \ ^ I 


Wtp (Wfc, Si, P/j + 0 ( ^ Sz - ^ (sfy - 2 s\"\si - sf^) 


1=1 


1=1 


1=1 


eev, 


(31) 


whieh leads to an upper bound of (1291) . Then, to tighten the obtained upper bound, we employ 
an iterative algorithm whieh is summarized in Table HIl First, we initialize the value of sf'^ for 
iteration index i = 0. Then, in eaeh iteration, we solve (l3TI) for given values of s[*\ ef. line 3, 
and update with the intermediate solution s^, of. line 4. The proposed iterative method 
generates a sequence of feasible solutions with respect to (|29|) by solving the convex 

upper bound problem (l3T]) successively. As shown in lf3^ . the proposed suboptimal iterative 
algorithm converges to a local optimal solutiorO of (l29l) with polynomial time computational 
complexity. In fact, the proposed suboptimal algorithm benefits from the convexity of (l3TI) 
and different numerical methods can be used to efficiently solve (l3T]) . In particular, when the 
primal-dual path-following interior-point method is used with a proper choice of kernel(/barrier) 
function, cf. ifSTll . If38]l . the computational complexity of the proposed suboptimal algorithm 
is 0(L ^^^ (NtL)'^ ln((A^TT)^/e)) with respect to N^L for a given solution accuracy e > 0 
lf39]l . where O(-) stands for the big-0 notation. The computational complexity is significantly 
reduced compared to the computational complexity of an exhaustive search which is given by 
0{2^'^^{NT:Lln{NT^L/e)) with respect to Nt^L, i.e., cf. Figure [3l 
Remark 4: The proposed algorithm requires to be a feasible point for the initialization. 


i.e., for i = 0. This point can be easily obtained since the constraints in 

V. Simulation Results 


span a convex set. 


In this section, we evaluate the system performance of the proposed resource allocation designs 
via simulations. There are L = 3 FD radio BSs in the system, which are placed at the comer 
points of an equilateral triangle. The inter-site distance between any two FD radio BSs is 250 
meters. The uplink and downlink users are uniformly distributed inside a disc with radius 500 
meters centered at the centroid of the triangle. We set the constant weights for the downlink and 
uplink power consumption as p = Q = 1, Vj G Ku}. The penalty term f for the proposed 


®By following a similar approach as in the proof of Theorem 1, it can be shown that Rank(Wfc) 
adopted SDP relaxation. 


1 holds despite the 




TABLE III 


System parameters 
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Carrier center frequency and path loss exponent 

1.9 GHz and 3.6 

Multipath fading distribution and total noise variance, cr^ 

Rayleigh fading and —62 dBm 

Minimum required SINK for uplink user j, 

10 dB 

Power amplifier power efficiency and antenna power consumption in idle mode, 

I/sd = 1/eu = 0.2 and 0 dBm 

Max. transmit power for downlink and uplink, and 

48 dBm and 23 dBm 


suboptimal algorithm is set to Also, Pq = 0 is adopted in all simulation result Unless 

speeified otherwise, we assume 50 dB of self-interference cancellatioiu at the FD radio BSs and 
the circuit power consumption per antenna is = 30 dBm. The antenna gains for the BSs 

and the users are 10 dBi and 0 dBi, respectively, and there are A^t = 20 antennas equipped 
in each FD BS resulting in = 60 antennas in the network. Furthermore, all downlink 
users require identical minimum SINRs, i.e., F^^^ = F^q, V/c. The performance of the proposed 
algorithms is compared with the performances of the following four baseline systems designed 
for peak system load when all the available antennas are activated. In particular, we minimize 
the total system power consumption of all four baseline systems using a similar approach as for 
the schemes proposed in this paper but set s; = 1, V/ G {1,..., N-^L}. The baseline systems are 
configured as follows. Baseline 1: a FD distributed antenna system (FD-DAS); Baseline 2: a HD 
distributed antenna system (HD-DAS); Baseline 3: a FD system with co-located antennas (FD- 
CAS); Baseline 4: a HD system with co-located antennas (HD-CAS). For the HD communication 
systems, we adopt static time division duplex such that uplink and downlink communication 
occur in non-overlapping equal-length time intervals. In other words, both self-interference and 
the uplink-to-downlink co-channel interference are avoided. For a fair performance comparison 
between HD and FD systems, we set log 2 (l-l-rJ^q,) = 1/2 log 2 (l-l-r^q“™) and log 2 (l-l-r°q,) = 


l/21og2(l + r 


DL-HD'i 
req- ) 


such that the minimum required SINRs for the uplink users, F 


UL-HD 
req. 


, and 


downlink users, become Fj^g”™ = (1 -t- Fj^q^^ — 1 and r°q“™ = (l-fF^q)^ — 1, 

respectively, to account for the penalty due to the loss in spectral efficiency of the HD protocol. 
Also, the power consumption of downlink and uplink transmission in the objective function of 
the HD systems is reduced by a factor of two as at a given time either uplink or downlink 
transmission is performed. For the CAS, we assume that there is only one BS located at the 
center of the system, which is equipped with the same number of antennas as all FD BSs in the 


’We note that the value of Pq does not affect the resource allocation algorithm design. 

*We assume a balun analog circuit is implemented in the FD radio BSs which can cancel 50 dB of self-interference CD. 
The residual self-interference is handled by the beamforming matrix Wfc via the proposed optimization framework. 










Objective value of Average total system 

equation (29) (dBm) power consumption (dBm) 



Fig. 2. Convergence of the proposed iterative algorithms. Fig. 3. Computational complexity versus the total number 

of transmit antennas in the system, NtL. 


distributed setting combined, i.e., Furthermore, for all baseline systems, we remove the 

maximum transmit power constraints imposed for the downlink and uplink transmissions, i.e., 
constraints C3 and C4. The key parameters adopted in the simulations are provided in Table Uni 
A. Convergence and Computational Complexity of the Proposed Iterative Algorithms 


Figure [2] illustrates the convergence of the proposed optimal and suboptimal algorithms for 
different minimum required SINRs for downlink users, F^^. There are = 4 downlink users 
and K\j = 2 uplink users in the system. It can be seen from the upper half of Figure |2] that 
the proposed optimal algorithm in Table I converges to the optimal solution in less than 350 
iterations, i.e., the upper bound value meets the lower bound value. On the other hand, from the 
lower half of Figure [2l we observe that the suboptimal algorithm converges to a local optimal 
value after less than 20 iterations. In the sequel, we show the performance of the suboptimal 
iterative algorithm for 20 iterations. 

Figure [3] compares the computational complexity of the brute force approach with that of the 
proposed suboptimal algorithnj^ for 20 iterations and solution accuracy e = 0.1. The system 
setting is identical to the scenario in Figure [2] and the results are computed based on the big- 
O complexity analysis in Section IV. As can be observed, the proposed suboptimal resource 
allocation algorithm requires a significantly lower computational complexity compared to the 
brute force approach, especially for large numbers of antennas. 


^ The proposed optimal algorithm may have the same computational complexity as the brute force approach in the worst case 
scenario although this seldom happens in practice. 






























Fig. 4. Average total system power consumption (dBm) versus the minimum required SINRs for the downlink users, r°q, for 
different systems. The double-sided arrows indicate the performance gain of the proposed FD system compared to traditional 
HD communication systems. 


B. Average Total System Power Consumption 

In Figure IH we study the average total system power eonsumption versus the minimum 
required SINRs of the downlink users, F^q. There are K-q = 4 downlink users and iTu = 2 
uplink users in the system. It ean be observed that the average total system power eonsumption 
inereases gradually with F^q. In faet, as the QoS requirements of the downlink users beeome 
more stringent, a higher downlink transmit power is needed to fulfill the requirement. At the 
same time, the self-interferenee power inereases with the downlink transmit power. Thus, the 
FD radio BSs have to utilize more degrees of freedom for self-interferenee suppression, and 
as a eonsequenee, less degrees of freedom are available for redueing the total system power 
eonsumption. On the other hand, the proposed suboptimal iterative resouree alloeation algorithm 
offers praetieally the same performanee as the optimal algorithm for the eonsidered seenario. As 
ean be observed, the two proposed algorithms faeilitate signifieant power savings eompared to all 
baseline system arehiteetures (whieh aetivate always all available antennas), espeeially for low to 
moderate system loads, i.e., F^q < 21 dB. Indeed, aetivating all antennas may not be benefieial 
for the total system power eonsumption when the load of the system is relatively small, sinee 
in this ease, the power eonsumption eaused by an extra antenna eireuit outweighs the power 
reduetion for information transmission offered by the extra aetivated antenna. Nevertheless, the 
performanee gap between the two proposed algorithms and baseline system 1 diminishes as 












Fig. 5. Average total system power consumption (dBm) versus the number of downlink users for F^q = 21 dB and different 
systems. The double-sided arrows indicate the performance gain of the proposed FD protocol compared to the FID protocol. 

the minimum required SINRs for the downlink users inerease. In partieular, the BSs are foreed 
to transmit with high power to satisfy the more stringent QoS requirements when the number 
of aetivated antennas is small. As a result, the two proposed algorithms have to aetivate more 
antennas, ef. also Figure [6l for improving the power effieieney of the system whieh yields 
a similar resouree alloeation as baseline system 1. Additionally, the two proposed algorithms 
outperform HD baseline systems 2 and 4 by a eonsiderable margin. As ean be seen, in the HD 
systems, an exeeedingly large system power eonsumption is required to meet the more stringent 
minimum required downlink SINRs to eompensate for the spectral efficiency loss inherent 
to the HD protocol. Furthermore, the distributed antennas deployed in the proposed systems 
provide spatial diversity across the network which shortens the distance between transmitters 
and receivers. This accounts for the power saving enabled by the two proposed algorithms 
compared to baseline CASs 3 and 4. 

Figure [5] depicts the average total system power consumption versus the number of downlink 
users for a minimum required downlink SINR of F^q = 21 dB. There are = 2 uplink users 
in the system. It is observed that the average total system power consumption increases with the 
number of downlink users. As more downlink users request communication services from the 
system, more QoS constraints are imposed on the optimization problem in (fT^ which reduces 
the size of the feasible solution set and thus results in a higher total system power consumption. 
In addition, the two proposed resource allocation algorithms outperform all baseline schemes 
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Fig. 6. Average number of activated antennas versus Fig. 7. Average number of activated antennas versus 

Freq. circuit power consumption per active antenna. 

due to the adopted optimization framework and the distributed antenna arehiteeture. 

C. Average Number of Activated Antennas 

In Figure [6l we study the average number of aetivated antennas versus the minimum required 
downlink SINR, for different numbers of downlink users. It ean be observed that the 
average number of aetivated antennas inereases with inereasing minimum required SINR for 
the downlink users. Although aetivating an extra antenna for signal transmission and reeeption 
eonsumes extra power in the eireuit, i.e., pActive_pidie > q ^ larger number of aetivated antennas 
inereases the degrees of freedom of the system whieh is benefieial if the QoS eonstraints are 
stringent. Speeifieally, with more antennas, the direetion of beamforming matrix ean be 
more aeeurately steered towards downlink user k whieh substantially reduees the neeessary 
downlink transmit power to aehieve a eertain QoS. Moreover, the redueed downlink transmit 
power also deereases the self-interferenee whieh in turn reduees the required uplink transmit 
power. In faet, for a small number of aetivated antennas, the FD radio BSs are required to 
transmit with exeeedingly high power if F^^ is large. As a result, the FD radio BSs prefer to 
aetivate more antennas to improve the power effieieney of information transmission, when the 
eost of aetivating extra antennas is less than the assoeiated potential transmit power saving. On 
the other hand, it ean be observed that the proposed sehemes aetivate more antennas when more 
downlink users are in the system. In faet, the downlink eo-ehannel interferenee inereases with 
the number of downlink users. Furthermore, the eo-ehannel interferenee eannot be suppressed 
by simply inereasing the downlink transmit power for all downlink users. Thus, extra spatial 
degrees of freedom are benefieial for deereasing the system power eonsumption. 

In Figure |71 we show the average number of aetivated antennas versus the eireuit power 
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consumption per active antenna, (dBm), for different minimum required SINRs for the 

downlink users. It is expeeted that the FD radio BSs prefer to aetivate more antennas when 
the eireuit power eonsumption per antenna is small or the SINK requirements of the downlink 
users are demanding, since in this ease, the power savings achieved by activating extra antennas 
surpasses the eorresponding eireuit power eonsumption. On the eontrary, when the eireuit power 
consumption per antenna is high, the FD radio BSs become more eonservative in aetivating 
antennas sinee using a large number of antennas may no longer be benefieial to the overall 
system power eonsumption. 

VI. Conclusions 

In this paper, we formulated the resource allocation algorithm design for power effieient 
distributed FD antenna networks as a mixed combinatorial and non-convex optimization problem, 
where the antenna eireuit power eonsumption and the QoS requirements of the uplink and 
downlink users were taken into account. Applying the generalized Bender’s deeomposition, we 
developed an optimal iterative resouree alloeation algorithm for solving the problem optimally. In 
addition, a polynomial time eomputational eomplexity suboptimal algorithm was also proposed 
to strike a balance between eomputational eomplexity and optimality. Simulation results showed 
that the proposed suboptimal iterative resource allocation algorithm approaches the optimal 
performanee in a small number of iterations. Furthermore, our results unveiled the substantial 
power savings enabled in FD radio distributed antennas networks by dynamieally switching off 
a subset of the available antennas; an exeeedingly large number of activated antennas may not 
be a eost effeetive solution for reducing the total system power eonsumption when the QoS 
requirements of the users are not stringent. 

Appendix 

A. Proof of Theorem [7] 

We start the proof by rewriting the Lagrangian function of the primal problem in (fTVl) in terms 
of the beamforming matrix W^: 

Ky> Kj, 

/:(©,$) = ^Tr(AfcWfc)-^Tr((Zfc + ^^)Wfc)+A (32) 

k=i k=i ^ '■safe 

Kii N'-pL N'pL N'-pL 

and Afc = + ^ pR, + ^ ^ (^D 

71 771,71 ^ ^ ^ 

j^k 1=1 m=l n=l 
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A denotes the eolleetion of variables that are independent of W^. For eonvenienee, the optimal 
primal and dual variables of the SDP relaxed version of (fTVl) are denoted by the eorresponding 
variables with an asterisk superseript. By exploiting the Karush-Kuhn-Tueker (KKT) optimality 
eonditions, we obtain the following equations: 


ZlhO, al>0,yk, 


(34) 


z * w * = o , 


7,* — A* 


a 






pDL 

reqfc 


(35) 

(36) 


where in (l3^ is obtained by substituting the optimal dual variables into (l3^ . From (l35l) . 
we know that the optimal beamforming matrix is a rank-one matrix when Rank(Z^) = 
N'j'L — 1. In partieular, is required to lie in the null spaee spanned by Z^ for ^ 0. 
As a result, by revealing the strueture of Z^, we ean study the rank of beamforming matrix 
W^. In the following, we first show by eontradietion that A^ is a positive definite matrix with 
probability one. To this end, we foeus on the dual problem in (l22l) . For a given set of optimal 
dual variables, $*, and a subset of optimal primal variables, {s*, P^*, W™’”*, P^m,n^ the 
dual problem in (l22l) ean be written as 

minimize (37) 

WfeGH^T V / 


Suppose A^ is negative semi-definite, i.e., A^ ^ 0, then we ean eonstruet a beamforming 
matrix = rw^w^ as one of the solutions of (iJTl) . where r > 0 is a sealing parameter and 
Wfc is the eigenveetor eorresponding to one of the non-positive eigenvalues of A^. We substitute 
Wfc = rwfcW^ into (1371) whieh yields 

^ T T 

£(©,$) = 52(*»wf (ZJ + + A, (38) 

k=i k=i ^ '■eqfc 

'-V-^ 

<0 

Besides, eonstraint Cl is satisfied with equality for the optimal solution and thus ak > 0. 
Furthermore, sinee the ehannel veetors of the downlink users, i.e., hn^., Vfc G {1,..., K-q], are 
assumed to be statistieally independent, we obtain — r Y^k=i 3^^ {^k^k {'^k + ^ -CX) 

when we set r —)■ CX). Thus, the dual optimal value beeomes unbounded from below. Yet, the 
optimal value of the primal problem in (fTTI) is non-negative for T^q^ > 0 whieh leads to 
a eontradietion as strong duality does not hold. Therefore, for the optimal solution, A^ is a 
positive definite matrix with probability one and Rank(A^) = N-^L, i.e., A^ has full rank. 
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Then, by exploiting (|3^ and basie rank inequality results, we have the following implication: 

Rank(Z*) + Rank > Rank (Z* + = Rank(A*) = N^L 

^ Rank(Z^) > Nt^L — 1. (39) 

Furthermore, 0 is required to satisfy Cl for > 0. Thus, Rank(Z^) = — 1 and 

Rank(W^) = 1 hold with probability one. ■ 

B. Proof of Proposition [7] 

We start the proof by studying the solution of the SDP relaxed version of (fT71) via its dual 
problem in (l22l) . For a given set of optimal dual variables ^*( 7 ), we have 0(i) 

= argndn Wtp (^W^, s,, +fi{&,^{i)) 

K^i N'^L N'^L 

Pi Tr(WfcR,) + ^m,nQm,n “t“ ^m^nQm^n ^m^nQ.m^ni (40) 

k=l /=! m=l n=l 

where the first equality is due to the KKT conditions of the SDP relaxed version of (fTTI) . On 
the other hand, we can rewrite function G {1,...,i}, in ^25\i as si^k) 


Kj^ Nt^L N^L 

minmize UTp(Wk, si, pfj +/i(0, $( 7 )) + ^^ Tr(WfcR;) + y^ '^i^ra,nqm,r. 

k=l (=1 m=l n=l 

"I Nt:L Nt:L 

P^m,nqm,n ^m,nQm,n Z T ^ ^ ^ ^ ^m,n(^^n T 1) ^m,n^m ^m,n^n 


m=ln=l 


Ku N-yL 


E E ^'-p. 

k=l l=l 


DL 

max; 


(41) 


The difference between (l40l) and (|4T]) is a constant offset. Thus, 0(f) is also the solution for the 
minimization in the master problem in (l4T]) for the f-th constraint in (I24cl) . The same approach 
can be adopted to prove that the solution of (1231) is also the solution of (|2^ . ■ 


C. Proof of Theorem |2] 


We start the proof of Theorem [2] by using the abstract Lagrangian duality GTll . BOll . iHTll . In 
particular, the optimization problem in (|2^ can be written as 


minimize maximize £(0,0) 
©ex' <;6>0 


( 42 ) 
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where 

N^L Nt^L 

Z(0,0) = Wtp (w,, 5,, ^ 0 ( ^ 5, - ^ sf) (43) 

i=i 1=1 

and the dual problem of (|2^ is given by 

maximize minimize £(©,0). (44) 

4>>0 ©GX> V ’-ry 

For notational simplicity, we define 

f2(0) = minimize £(0,0). (45) 

Then, we have the following inequalities: 

maximize fl(0) = maximize minimize £(©,0) (46a) 

<p>o <p>o &ev 

(a) _ 

< minimize maximize £(0. 6) = (l2^ . (46b) 

©GX> 0>O 

where (a) is due to the weak duality Il3^ . We note that ~ J2^=i^ > 0 for © G P such 

that £(©, 0) is a monotonically increasing function in 0. In other words, fl(0) is increasing in 
0 and is bounded from above by the optimal value of (l4^ . Suppose the optimal solution for 
(146a| ) is denoted as 0 o and ©* = {W^, where 0 < 0 o < oo. Then, 

we study the following two cases for the solution structure of (I46al) . In the first case, we assume 
= 0 for (I46al) . As a result, ©* is also a feasible solution to (1281) . Subsequently, 
we substitute ©* into the optimization problem in (|2^ which yields: 

fi(0*) = Wtp (Wfc, Si, pu) > m- (47) 

By utilizing (I46al) and (1471) . we can conclude that 

minimize maximize £(©, 0) = maximize minimize £(©,0) (48) 

©GX> </.>0 <t>>0 &GV 

must hold for ~ '5/^ = 0- Furthermore, the monotonicity of 17(0) with respect to 

0 implies that 

17(0) = (HH), V0>00 (49) 

and the result of Theorem [2] follows immediately. 

Now, we study the case of ~ > 0 at the optimal solution for (I46al) . The 

optimization problem maximize 17(0) —)■ cxo is unbounded from above due to the monotonicity 
of function 17(0) with respect to 0. This contradicts the inequality in (I46al) as (l28l) is finite and 
positive. Thus, for the optimal solution, Si — -s? = 0 holds and the result of Theorem 

|2] follows immediately from the first considered case. ■ 
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